Scale Invariance and Self-averaging in disordered systems. 
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ABSTRACT 

In a previous paper [Q we found that in the random field Ising model at zero 
temperature in three dimensions the correlation length is not self- averaging near 
the critical point and that the violation of self-averaging is maximal. This is due to 
the formation of bound states in the underlying field theory. We present a similar 
study for the case of disordered Potts and Ising ferromagnets in two dimensions 
near the critical temperature. In the random Potts model the correlation length is 
not self-averaging near the critical temperature but the violation of self-averaging is 
weaker than in the random field case. In the random Ising model we find still weaker 
violations of self-averaging and we cannot rule out the possibility of the restoration 
of self-averaging in the infinite volume limit. 
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One of the important issues in the physics of disordered systems is the understanding of 
which physical quantities are self averaging, i.e. for which observables the sample to sample 
fluctuations vanish in the thermodynamic limit [IJ |3] ■ 

In a recent workup we studied the behavior of the correlation length (CL) of the random 
field Ising model (RFIM) in three dimensions at zero temperature. We found surprisingly that 
the CL is not self-averaging and that the violation of self-averaging is maximal in the vicinity 
of the critical point. This is due to the formation of bound states in the underlying replica 
field theory. The interaction among replicas is attractive and therefore it has the capability 
of forming bound states. In the formation of bound states there is competition between the 
strength of the attractive forces and the size of the available phase space. We have argued that 
in two dimensions phase space is small and even a small attraction should win. 

In this paper we address the question of the existence of bound states and their consequence, 
the non self-averaging of the CL, in the case of two disordered ferromagnetic spin models in 
two dimensions: the q = 3 states Potts model and the Ising model. It is known that, in 
agreement with the Harris criterion[2j, the disorder is relevant for the q = 3 states Potts model, 
i.e. the renormalization group flows to a stable fixed point, which is different from the fixed 
point of the pure system^ Ej. For the Ising model the disorder is marginally irrelevant and 
the renormalization group flows back to the fixed point of the pure system with logarithmic 
corrections (ZJ. 

In the case of the Potts model we find very similar behavior as in the RFIM. A scaling 
analysis shows that, as in the RFIM, the violation of self-averaging of the CL, although weaker 
than in the RFIM, persists for arbitrarily large systems in the vicinity of the critical point. 

In the case of the Ising model, the violation of self-averaging is still weaker. The scaling 
analysis seems to work but we cannot exclude a very slow restoration of self-averaging in the 
thermodynamic limit (see later). 

We conclude that the non self-averaging of the CL is not a singularity of the RFIM at 
zero temperature, but a more general phenomenon in the physics of disordered systems. The 
correlation length plays a crucial role in the scaling theories of phase transitions. It is not 
clear yet what are the consequences of the absence of self- averaging of the CL for those scaling 
theories in disordered systems. 

In both the Potts and the Ising models we consider the case of square lattices of size I? 
where L = 200, 350, and 500, with nearest neighbor interactions. The ferromagnetic couplings 
J are independent random variables which can take the two values J = J\ and J = J 2 with 
equal probability. It is well known that these models are self-dual and that the critical inverse 
temperature j3 = j3 c is given by |Hj: 

(exp/3 c J 1 -l)(exp/? c J 2 -l)= q. (1) 

The knowledge of (5 C greatly simplifies the scaling analysis. It is also known that cluster 
algorithms|2] are very effective for these models. We simulated 1000 samples for L = 500, 
2000 for L = 350 and 10000 for L = 200. As in the case of the RFIM[T], we measured the 
correlation length by studying the dependence on the boundary conditions. More precisely in 
the Ising case (the generalization to the Potts case is obvious) we set all the spins equal to one 
on the line x = and choose free boundary conditions at the other end of the lattice x = L. 
We impose periodic boundary conditions on the perpendicular direction. After thermalisation, 
for every sample s we measure the magnetization m s (x) = l/LJ2 y < &{x,y) > and study its 
dependence on x and on the inverse temperature (3. For the Potts model we studied extensively 
the case Ji/J 2 = 10jm and P/Pc = -95, .96, .97, .98, .99, .995 and 1. We verify thermal- 
isation in the following ways. For every realization of the couplings, we simulate two copies 
with different initial conditions. In one copy all spins are set to +1 and in the other they are 
chosen at random (except on the line x = 1, where they are set to +1). Thermalisation time tth 



1 



0.05 



c/5 
05 



0.045 - 

0.04 - 
0.035 - 

0.03 
0.025 

0.02 - 
0.015 - 

0.01 - 
0.005 - 




m 

_i_ 



0.95 0.96 0.97 0.98 0.99 1 

Figure 1: Mass as a function of f3 < j3 c for the q = 3 Potts model and L = 500. 



must be long enough in order for the magnetization m s (x) of the 2 copies to be indistinguish- 
able. The local magnetizations < cr(x, y) > are the time averages of the local spins during the 
measuring time t m We measure the autocorrelation time r and require both the thermalisation 
time tth and the measuring time t m to be much larger than r. We also verify that there are 
no time drifts in the measurements, i.e. averages over the first half of t m do not differ from 
measurements over the second half of t m . Even at /3 = (3 C and for the volumes we simulated, 
tth is finite and reasonably small because of the cluster algorithm. For q = 3, we measured 
tl=200 — 200, t l=350 ~ 400, t l=50 o ^ 800 and we choose t th = 10r, t m = 20r. 

To reduce errors, we found convenient to use a variance reduction technique: we measured 
the probability p(x,y) of the spin a(x,y) to be on the same cluster with the spins at the 
boundary x = 0. Obviously < a(x,y) >= 2p(x,y) — 1. In this way we eliminate statistical 
fluctuations due to flips of large clusters. 

In order to study sample to sample fluctuations, we measured, as in the case of the RFIMpQ, 
m(x) = m s (x) and m( 2 )(x) = m s (x) 2 , where, as usually, the bar denotes average over the 
random coupling samples. As it was mentioned earlier, for every sample we simulated two copies, 
differing only at the initial conditions. In order to reduce errors, we measure the probability 
P2(x,y) that in both of the two copies we simulated, the spin at the site (x,y) belongs to 
the cluster connected to the boundary at x = 0. Obviously < ai(x,y) >< d2(x,y) >= 
2p2(x, y) — 1 where a\ and a"2 are the spins of the two copies. With these definitions, m(o\{x) = 

For j3 < (3 C , we expect m(x) and m^ 2 )(x) to be exponentially decreasing functions of x. We 
fitted m(x) with the formula 



exp (-fix) b 

mix) = (oH — 

v 1 x a v x' 



(2) 



We found that this formula fits very well the data for all values of x except near the boundaries 
x ~ and x ~ L and this for all values of (3 we studied. The mass fx (or inverse correlation 
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Figure 2: Finite size scaling plot of the mass ratio (see text) versus L 1/,!/ (/3 c — /3)//3 c for the 
q = 3 Potts model. 



length) for the q = 3 states Potts model, J\j Ji = 10 and L = 500 is plotted in figure ^ as a 
function of (3 < (3 C . Throughout this paper we used the jack-knife method to estimate statistical 
errors. We found that statistical errors are smaller than systematic errors. We applied a small x 
cutoff xq and a large x cutoff x\ in fitting the data, because, as expected, close to the boundaries 
the previous exponential form is not valid. The choice of xo and x\ is an important source of 
systematic errors on the mass fi. For every size L and every /? we tried all pairs of xo and x\ 
for which the fit was reasonable, i.e. the x 2 °f the fit per degree of freedom was .10 < x 2 < -80 
(the values of m(x) at different x are not statistically independent). The error bars in figure ^ 
reflect the dependence of the mass fi on the different choices of "acceptable" cutoff's (statistical 
errors are much smaller). Typical values were 5 < xo < 15 and .5L < x\ < .8L depending on 
L and (3. We performed a similar analysis of m(o\ (x) defined above and measured in the same 
way the corresponding mass The ratio R = ^2/^ is dimensionless. It is plotted in figure 
121 as a function of the scaling variable L 1//l/ (/3 c — /3)//3 c for the q = 3 states Potts model and 
J1/J2 = 10. Following the results of previous simulations ^H], we assumed that the critical 
exponent v = 1. 

Figure [21 shows that the behavior of R is compatible with the scaling ansatz. For (3 far from 
P c R = 2, i.e. for large x, m(2)(x) decays twice as fast as m(x). This shows that, far from 
criticality, the correlation length is self-averaging. 

Near the critical point, R < 2. As this is discussed in this implies the existence of a 
bound state in the underlying field theory. In the case of the RFIM R = 1, which means that 
the bound state completely dominates the asymptotic long distance behavior. In the q = 3 
random Potts model 1 < R < 2, which suggests a weaker coupling of the bound state, as we 
will see. 

As this is discussed in pQj if there exist more than one comparable mass scales, as it can 
be in the case of bound states, we expect a superposition of exponentials governing the long 
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Figure 3: Finite size scaling plot of c\ (see text) for the q = 3 Potts model. 



distance behavior. This was the case for the RFIM. 

We performed a similar analysis and fitted our data as follows: 

m (2)( x ) = cim(x) + C2(m(x)) 2 . 



(3) 



It turned out that this two parameter fit of m,2(x) is better than the four parameter fit of 
eq. @. If c% 7^ the first term dominates for large x (m(x) decreases much slower than m(x) 2 ) 
and the mass is non self-averaging. As for the RFIM we made the finite size scaling hypothesis 
that the dimensionless coefficient c\ is a function of L x l v ((3 C -/3)//3 c . In figure|3]we plot ci versus 
HPc ~ 0)/Pc (we took v = 1) for different values of f3 and L = 200, 350 and 500. We see that 
for L = 500 and L = 350 the data fall completely on each other, while for L = 200 small non 
leading corrections to scaling seem to be present. For j3 = /3 C , c\ = .10 ± .02. We conclude that 
for the q = 3 random Potts model the mass is not self-averaging for any finite volume, provided 
the temperature is close enough to its critical value. This is due to the presence of bound states 
in the replica field theory, as in the case of the RFIM. The violation of self-averaging is much 
weaker than for the random field model in three dimensions where we found c\ = 1, i.e the 
maximum possible violation of self-averaging. 

We performed a similar analysis for the Ising model with J1/J2 = 5. In figure |I] we plot c\ 
versus L(f3 c — (3)/(3 c . At first sight the data seem again compatible with this scaling hypothesis 
and there is again evidence of violation of self-averaging of the mass, although weaker than in 
the Potts case ( c\ ~ .033 at f3 c ). A closer look to the data at j3 = (3 C , however, leaves open 
the possibility of a slow crossover to c\ = at (3 C as L — > 00. (Away from f3 c scaling seems to 
work fine.) In order to settle this question we simulated a larger lattice with L = 700 at (5 C . 
If scaling of c\ holds, c\ should be independent of L at /3 = f3 c . We found the following values 
of ci at f3 c : .034 ± .003 for L = 700, .0355 ± .004 for L = 500, .0365 ± .004 for L = 350 and 
.0425 ± .004 for L = 200. c\ seems to decrease from L = 200 to L = 350 and then change very 
little from L = 350 to L = 500 and to L = 700. It is known j7j that the renormalization group 
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Figure 4: Finite size scaling plot of c\ (see text) for the Ising model. 



fixed point of the random Ising model flows to the pure Ising fixed point at a logarithmic rate. 
Our data cannot discriminate between such a behavior or sub-dominant corrections to scaling 
when the lattice is not large enough. 

Finally in order to verify the dependency of the violations of self-averaging of the CL on the 
disorder, i.e. the ratio r = J1/J2, and the number of states of the Potts model q, we measured 
ci for L = 200 and (3 = (3 C for other values of q and r. For q = 8 and r = 10 we got c\ = .19± .02. 
For q = 2 (Ising) we got ci = .0145 ± .005 for r = 2 and ci = .063 ± .002 for r = 10. We 
conclude that c±, i.e. the violations of self-averaging of the correlation length, increase when 
the disorder is larger, or when the number q of states of the Potts model increases. 
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